APPLIED 

MATHEMATICS  INSTITUTE 


University  of  Delaware 
Newark,  Delaware 


R  9.  02 


8  :M2 


Approv'^.l  f  O’"  pnMlc  rolease 

diutrilution  ’onliniited. 


iUU' 


d) 


ON  MODIFIED  GREEN'S  FUNCTIONS 
IN  EXTERIOR  PROBLEMS 
FOR  THE  HELMHOLTZ  EQUATION 


R.  E.  Kleinman 

Department  of  Mathematical  Sciences 
University  of  Delaware 
Newark,  Delaware  19711 


G .  P .  Roach 

Department  of  Mathematics 
University  of  Strathclyde 
Glasgow  G1  IXH,  Scotland 


APPLIED  MATHEMATICS  INSTITUTE 
Technical  Report  #  116A 


December ,  1981 


air  F 


J.  • 
Tocii.-.  1, 


Vct-r- 

r  rue  TM- 

’  Just:;  ■  ,  ' 


(AFSC)  ! 

I 

■  ^nd  is  I - 

I  AV'  i  1 


*  on  1  i V  i3  ioa 


i 


I  Introduction 


Recently,  Jones  has  presented  a  method  for  overcoming  the 
non  uniqueness  problem  arising  in  boundary  integral  equation  formu¬ 
lations  of  the  Dirichlet  and  Neumann  problems  for  the  Helmholtz 
equation.  The  major  portion  of  Jones'  analysis  concerned  the  exterior 
Neumann  problem  in  two  dimensions  although  he  indicated  hov  the 
results  generalised  to  three  dimensions  and  suggested  that  the  exterior 
Dirichlet  problem  could  be  similarly  treated.  I’rsell  simplified 
the  proof  of  a  k.ey  theorem  in  but  confined  his  remarks  to  the 
exterior  Neumann  problem  in  two  dimensions.  In  .3_  the  authors 
presented  a  systematic  exposition  of  boundary  integral  equation 
formulations  of  both  Dirichlet  and  Neumann  problems  and  presented  a 
number  of  useful  properties  of  the  boundary  integral  operators  arising 
in  both  layer  theoretic  and  Green's  function  approaches. 

in  particular  it  was  shown  that  uniqueness  of  the  boundary  integral 
equation  formulations  of  exterior  problems  could  be  retained  even 
at  eigenvalues  of  the  corresponding  adjoint  interior  problems  by 
creating  a  pair  of  coupled  equations.  That  treatment  dealt  with  three 
dimensional  problems  although  the  results  remain  intact  when  the 
fundamental  solution  of  the  Helmholtz  equation  in  n  dimensions  is  used. 

In  the  present  note  we  show  how  Jones'  modification  can  be 
incorporated  into  the  boundary  integral  equation  formalism  of  _3_  . 
Ursell's  simplification  is  adapted  to  three  dimensions  and  explicit 
results  are  obtained  for  both  Dirichlet  and  Neumann  problems.  In 
particular  it  is  demonstrated  that  a  single  boundary  integral  equation 
is  uniquely  solvable  in  each  case  even  at  interior  eigenvalues  of  the 
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adjoint  problems  by  suitably  modifying  the  Green's  function  in  a  way 
suggested  by  Jones'  approach.  Moreover  in  the  present  note  it  is  shown 
that  by  abandoning  the  restriction  to  real  coefficients  in  the 
modification  which  Jones  and  TJrsell  found  sufficient  to  eliminate  non 
uniqueness  at  interior  eigenvalues,  the  coefficients  may  be  chosen  to 
be  optimal  with  respect  to  certain  specific  criteria.  In  particular 
results  are  presented  which  show  how  to  choose  the  coefficients  so  as 
to  minimize  the  difference  between  the  modified  and  exact  Green's 
functions  for  the  Dirichlet  and  Neumann  problems  and  furthermore  an 
algorithm  is  presented  which  determines  the  coefficients  so  as  to 
minimize  the  spectral  radius  of  the  modified  boundary  integral 
operator.  Different  coefficient  choices  result  in  each  case. 

An  alternate  treatment  involving  the  explicit  determination  of 
a  set  of  eigenfunctions  of  the  unmodified  operator  is  given  in  . 

II  Notation  and  definitions 

We  adopt  the  notation  of  [3}  for  present  purposes.  Thus  let 
D_  denote  a  bounded  domain  in  ]R  *  with  boundary  D  and  exterior  D_^. 

The  boundary  3D  will  be  assumed  Lyapunov  and  when  more  stringent 
smoothness  assumptions  are  needed  they  will  be  explicitly  noted. 

Let  R  =  R(P,Q)  denote  the  distance  between  any  two  points  P 
and  Q  in  ]R  - .  A  function  u^(P)  which  satisfies 

(V^  +  k2)u_^(P)  -  0,  P  £  D_^U  3D,  (1) 

3u 

lim  r  ( - iku  )  »  0 

r  P  5r 
P  P 


will  be  called  (following  Wilcox  j_5' 
Helmholtz  equation  where 


)  a  radiating  solution  of  the 
are  the  spherical  polar  coordinates 
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of  a  point  P  relative  to  a  Cartesian  coordinate  system  with  origin  in 
D_.  It  should  be  stressed  that  the  notation  in  (1)  is  meant  to  imply 
that  there  is  an  open  domain  including  the  boundary  3D  in  which  is 
differentiable.  A  fundamental  solution  of  the  Helmholtz  equation  is 
given  by 

ikR 

y„(P,Q):  -  -  ■= -  (5) 

2:tR 

and  if  g(P,Q)  is  a  radiating  solution  of  the  Helmholtz  equation  in 
both  P  and  Q  then 

>^(P,Q):  -  YqCP.Q)  g(P.Q) 

_s  also  a  fundamental  solution.  If  w  t  L2(3D)  then  we  may  obtain 
standard  and  modified  forms  of  single  and  double  layer  distributions 
of  density  w  according  to  whether  or  is  employed  as  the 
fundamental  solution.  Thus  we  have,  for  j  *•  0,1, 


(SjW)(P): 


w(q)  Y. (P,q)ds  ,  P  e  IR  - 
5D  ^  ^ 


(D.w)(P):  -  ;  w(q)^  (P,q)ds  ,  P  E  IR-  (6) 

J  -3D 

where  is  the  derivative  in  the  direction  of  the  outward  normal  tc 
3n 

q 

3D  at  the  point  a  (n  points  into  D  )  .  Further  we  write  —  and 
*  q  +  on  + 

5  P 

to  denote  the  normal  derivatives  when  P  -  p  e  5D  from  D  and  D 
on~  +  - 

P 

respectively. 

Denote  by  Kj  the  boundary  integral  operator,  compact  on  L.,(3D)  as 
well  as  C^(3D),  (Mikhlin  [bj ) ,  defined  by 


(K.w)(P):  =  :  w(q)  (P,q)  ds  ,  p  e  5D 

J  -  3D 


with  L,  adjoint 


*  f  ‘’^1 

CK.  w)(P);  --  1  w(q)  (P,q)  ds  ,  p  £  3D  (8) 

J  '  3D  ^""q  ^ 

where  a  bar  is  used  to  denote  complex  conjugate.  In  terms  of  the 
jump  conditions  for  the  single  and  double  layer  distributions  are 


(S^w)  •=  +  w  +  K^w,  p  E  3D  (9) 

3n” 

P 


and 


lim  D.w  =  +  w  + 
P-^p+ 


K. 

J 


w, 


p  £  3D , 


(10) 


Note  that  g(?  Q)  must  be  defined  in  a  neighbourhood  of  3D  in  rder 
for  (9)  to  have  meaning. 

If  the  normal  derivative  of  the  double  layer  exists  (e.g.  if 
w  £  Cj^OD))  then,  e.g.  Gunter  [7j  , 


(DjW) 


(11) 


The  existence  of  the  normal  derivative  of  the  double  layer  when 
the  density  arises  as  the  solution  of  an  integral  equation, as  in  the 
present  context,  is  established  by  Angell  and  Kleinman  [SJ  .  Krai  [9_ 
discusses  the  case  of  even  weaker  hypothesis  on  the  density. 

If  u_^  is  a  radiating  solution  of  the  Helmholtz  equation  then 
Green's  Theorem  tog-ether  with  (9)  -  (11)  yield  the  representation 


/  3u  \ 

-_)(P)  -D^u^(P) 

\  cn  / 


HZu^  ,  P  t 

j  ,  ?  E  3D 

0  ,  P  E  D_ 


(12) 


Since  g(?,q)  must  have  singularities  in  D_  (there  are  no  nontrivial 
radiating  wave  functions  in  the  representation  (12)  does  not 
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always  hold  if  is  replaced  by  However  the  following  is  valid: 

/3u, 


2u^(P),PcD 
\3n  /  >  p  c  aD 


(13) 


The  representations  (1?)  and  13)  yield  the  boundary’  integral  equations, 


-* 


K 


I - r!  -  T-  (D.  u.)  *  - 


^  \3n^i  -  3n 


p  £  3D. 


(U) 


(13) 


Observe  that  y^  is  an  approximate  Green's  function  as  defined  by  Roach 

1“  “-I 

J.0,  11_,  and  the  present  analysis  may  be  considered  an  application  of 
the  ideas  presented  there. 

Finally  we  adapt  Ursell's  notation  to  IP  •  and  define  normalized 
spherical  wave  functions 


& 

C  ’*(P); 

nm 


(2n+l)  Z®’^  (kr  )p“(cos  3  )cos  me 

In+m)  .J  n  P  n  P  P 


r  ik 

L  2r  Sn 

S  ®'i(P):  c  (2n.l)  V«^(kr  )p“ 

nm  L  2ir  m  (n+m)TJ  n  p  r. 

fl,  m  =  0 

"  « 


(16) 


and 


where 


(cos  f  )  sin  me 
P  P 


2,  m  >  0 

V/ 

P™  (cos  6  )  is  the  associated  Legendre  pol^/nomial  and  (kr)  are 

n  p  "  n 

the  spherical  Bessel  and  Hankel  functions 


Z  (kr)  *  h^*  (kr) 
n  n 


(kr)  =  (hr). 


Similarly  we  denote  differentiation  with  respect  to  kr^  by  a  prime  so 
chat 


:e,i' 

(P): 

3 

(?)  . 
nm 

nm 

:e,i' 

t 

(?) : 

S”’  (?)  , 

'nm 

\  *  /  • 

3(krp) 

nm  '  ’ 

(17) 
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and  che  differentiation  only  affects  the  radial  functions. 

The  orthogonality  of  spherical  harmonics  leads  to  the  relations 


3B 


C  (P)S®’^(P)ds 
nm  '  JLs 


^  ]  C®’^(P)S®’^’cP)ds^  -  0  yn.m.i.s.  (18) 

P  j  jB  P 

’^P 


3B 


c®;\p)c!!'-  (p)ds^ 

ntn  jk.S  p 


0,^  #  n  and/or  n  #  s 


2  k  r^  Z^^’Ckr^)  2®’*  (hr^) , 

p  n  p  n  p 


n  ”  A, ,  10“  s  < 


(19) 


and 


SB. 


S^^(P)S®’^  (P)ds  “  0,x  f  n  and/or  m  #  s  or  m  =  0 
nm  is  p 


2  k  r2  Z®’^(kr^)  2®’^'  (kr  ), 
p  n  P  n  ? 


n  =  m  ®  s . 


(20) 


where  3B  is  Che  surface  of  a  sphere  of  radius  r  =  'p  . 

^P  ^ 

Because  of  the  expansion  of  the  fundamental  solution  which 

notation  is 

*  n  . 

Y  (p,Q)  -  :  -  {c^(p<)c®„(p>)  +  s^^(p<)s®  (?>):■ 

o  ^  nti  nm  nm  nm 

n“0  n«0 


in  tnis 


wnere 


P  <  .  fp  i£  rj,  <  rn  P  >  -  ?  if  rj,  > 


Qifr,  <rp  |  .  jQi 


’■q  "  '^P. 


the  single  and  double  layers  have  che  representations 
®  n  . 

(S  w)(P)  -  ;  :  C®  (P)  +  S®  (?) ,  r  <  min  r 

o  n  n  nm  nm  P  ■vt,  <1 

n*0  m“0  qt3D 


(21) 


m  n 


;  :  C"  C®  (F)  *  S®  (P),  r  ^  max  r 

-0  m-0  ^  qc5D 


(22) 


and 


Illl  IIKJlIliMWIil.i 


T 


(D^w) (P) 


n*0  iD*0 
n 


I  C®  (P)  +  S®  (P) ,  r  4  min  r 

‘^nin  tun  ^nin  nxn 


^  qcoD 


*  f  *  I 

I  :  C®„  (P)  +  S®  (P)  ,  r  max  r 


where 


n»o  m»0 


C®’^ 

^nm 


C®’^ 

^om 

S®’^ 

oiiin 


S®’^: 


r^xm  ntn 


njD 


^  qcoD 


(23) 


''3D  -  q 

III  The  Exterior  Neumann  Problem 


:  w(q)  (q)  ds  , 

J  3D  ^ 

'  3D  ^^q  ™  ° 

[  q(q)  (q)  ds^. 

'  aiJ 


!  ^  ‘^"q 


(24) 


The  exterior  Neumann  problem  consists  of  finding  u_^,  a  radiating 
solution  of  the  Helmholtz  equation  in  the  sense  of  (1)  anc.  (2)  for 


which  =  g_^,  p  £  3D  where  g_^  is  a  given  function  on  the  boundary. 
Boundary  integral  equation  formulations  are  obtained  from  Green's 
theorem  (14)  and  (15). 

(I  +  K>)  w*  =  Sj  g^  (26) 

“*  •  «■) 

in  which  case  the  solution  is  represented  in  D_^  via  (12)  and  (13)  as 

(28) 


1 


%  *  f  -  2  P  -  D+- 

Alternatively  the  layer  theoretic  approach  starts  with  the  assumed 
form  of  the  solution 

u_^  »  (S^w)  (P) ,  P  £  D_^u'  3D 

which  with  (9)  yields  the  boundary  integral  equation 


(29) 


•  * 


8. 

(I  +  Kj)  w  -  g^.  p  e  3D  .  (30) 

As  in  [3!  we  denote  chose  real  values  of  k  for  which  the  homogeneous 
equations 


(I  +  Kj)  w 


(31) 


has  nontrivial  solutions  as  characteristic  values  of  (-K.)  and  obser'>/e, 

J 

since  Kj  is  compact  on  L2(3D),  that  if  k  is  a  characteristic  value  of 
(-Kj)  it  is  also  a  characteristic  value  for  (~K|) ,  (-K. )  and  (-Kt) . 

It  is  of  course  these  characteristic  values  which  impede  the  solution 
of  (26)  and  (30) .  In  [3j  it  was  shown  that  the  pair  of  equations 
(26)  and  (27)  had  a  unique  solution  for  all  real  values  of  k  when 
j  =  0.  The  fact  chat  a  pair  must  be  considered  coupled  with  the 
complicated  nature  cf  (27)  makes  this  approach  difficult  to  follow  in 
actually  constructing  solutions,  for  example  by  numerical  methods. 

It  is  this  complication  that  is  avoided  by  properly  perturbing  the 
fundamental  solution. 

First  we  note  the  following  important  result. 

Theorem  3.1:  For  Im  k>0 

(I  w  =  0  if  and  only  if  Sj  w  =  0,  pc  5D. 

Proof:  If  w  *  0  the  theorem  holds  vacuously.  If  j  =  0  the  result  is 

proven  in  ' .  For  j  ®  1  the  same  proof  applies.  StW=0,  p  c  3D 

implies  S,w  =  0,  ?  c  from  the  uniqueness  of  solutions  of  the 

exterior  Dirichlet  problem  for  Im  k  5.  0.  Hence  — ^  S^w  =  0  which 

3n 

with  (9)  shows  that  (I  +  K^)  w  *  0.  Conversely  if  (I  +  K)  w  =  0  then 
again  Sj^w  =0,  P  c  D_^  t;  3D  otherwise  it  would  violate  the  uniqueness 
theorem  for  Che  exterior  Neumann  problem  which  completes  the  proof. 

Also  we  have  for  nontrivial  w 
Theorem  3.2:  For  p  c  oD, 
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(I  +  K^)  w  -  0  <■“>  ■■  0  <•“>  k  is  an  eigenvalue  of  the  interior 

Dirichlet  problem. 

Actually  -k^  is  the  eigenvalue  of  the  Laplacian  but  in  the  present  paper, 
as  in  [3],  we  shall  understand  by  eigenvalues  of  the  interior  Dirichlet 

(Neumann)  problem  those  values  of  k  for  which  there  are  non  trivial 

•  ^  du 

solutions  of  +  k^)  u  *  0,  P  e  D  and  u  *  0  (•r—  =  0) ,  p  c  3D. 

“  3n 

This  result  is  proven  in  [3^  .  The  fact  that  characteristic  values  of 
-K^  are  identical  with  eigenvalues  of  the  interior  Dirichlet  problem 
does  not  generalize  to  because  g(P,Q)  is  not  defined  throughout 
D_.  Nevertheless,  with  a  suitable  choice  of  g(P,Q)  we  may  establish 
Theorem  3.3:  If 


g(P.Q)  = 


n“0  m”0 


{a  C®  (P)C®  (Q) 

nm  nm  nm 


b  „  S®  (P)S®  (Q)} 
nm  am  na 


either  a  “  0  or  .2a  ♦I  '<1  (33) 

nm  nm 

and  either  b  =0  or  :2b  +1  i<l  (34) 

nm  nm 

and  k  is  a  characteristic  value  of  -K^  then  k  is  an  eigenvalue  of  the 
interior  Dirichlet  problem. 

Proof:  Assume  that  k  is  a  characteristic  value  of  i-e* 

(I  +  K^)  w  -  0  (35) 

has  a  nontrivial  solution.  Then  Theorem  3.1  implies  that  S^w  *  0, 

P  t  D^^CjSD  and  since  g(P,q)  is  an  analytic  function  of  P  for  ’t  0, 
it  follows  that  Sj^w  »  0  is  defined  for  ?  t  D_s''0-  and  continuous  on  3D, 
Thus  defining 

u_:  -  S^w,  P  £  D_\{C}  (36) 


it  follows  that  u_  (hence  u_)  is  a  solution  of  the  Helmholtz  equation 
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in  D  n{0}  which  vanishes  on  3D.  Hence 


3u  3u 


(u  - —  ~  u  T — )  ds  *  0. 
3D  ■  ‘  P 


(37) 


Let  r  be  Che  radius  of  a  sphere  B  entirely  in  D  (r  $  min  r  ) . 
P  tp  p  '  q 


Then  with  Green's  theorem  the  equation  (37)  may  be  written 


r 

J  SB 


3u_  _  3u_  r  ,  -  — ~ 

(u  - - u  -r — )  ds  *  0  *  I  (S.w  S,w  -S.w  —  S  v)ds 

~  3n  -  3n  p  I  1  dn  1  J-  an  1 

r  P  P  p  P 

P  P 


(38) 


No  volume  integral  term  arises  because  k  is  a  characteristic  value  of 


-K,  ,  hence  real .  On  B 
■L  r 


S  w  =  j  g(P,q).w(q)  ds  = 
®  J  3D  ^ 


=  :  Z  {a  C®  C®  (P)  +  b  S®  S®  (P)}  (39) 

.  «  nm  ^nm  nm  nm  'x.nm  nm 

n*0  m“0 


where  we  use  the  notation  of  (24)  and  (25)  with  w  in  place  of  w. 
Substituting  (22)  and  (39)  in  (38)  we  find 


.  CO  n 

k  j  { _ 

'  3B  n^O  m*«0 
r 

P 


Z  (c®  (c^  (P)  +  a  C®  (P)}  +  S®  {s^  (P)  +  b  S®  (P) ;)  x 

V'uun  nm  nm  nm  'Uim  nm  nm  run  / 


jc»0  s*0 


C®  (C^  (P)  +  a„  C'  (P)}  +  S'  (S:  (P)  +  b.  S'  (?);)- 
'■uA.S  is  is  is  '\,is  JCS  A,S  JCS  ' 


;:e 


e  '  r. 


_  ;  T 


-p' 


00 

-  Z  Z  fc®  {c"‘  (P)  +  a  C"  (?)}  +  S'  -S^  (P)  +  b  S'  (P)})  x 

^  \.vnm  nm  nm  nm  ..^nm  nm  nm  nm 


:e  rri 


re  .ri 


re 


n*0  m*0 


i=0  s=0 


Z  (c®  {C^'(P)  +  a  C® 


®'(P)}+  S®„{S^'(P)  +  b,^  S®^  (P):]  }ds  -  0. 


s'  '  ^is  is'  is  is 


(40) 
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Upon  applying  (18)  -  (20)  this  becomes 


2(kr  )2  Z  z  {lC®„|2  (j  j'  +  a  j'  +  i  j  '  + 

P  jj_Q  -vnm'  '■-'n-^n  nm  n  -'ll  nm  ■'n  n 


I  2  >,(1)  u(2)\  ^ 

a  I  h  h  )  + 

'  nm'  n  n 

-  Is'  (j.j;  •  h'^)  j.  .  b  i  h<«'  .  ^b  :2h<»  - 

nm  n  n  nm  n  Jn  nm  "^n  n  nm  n  n 

-  !C'  |2  (j  ii  -  a  i  b«)'  *  ;  h(2>j.  .la  '^h<«b<’-)’)  - 

nm  n  n  nm  n  n  nm  n  n  nm  n  n  ^ 

■  ^  +  lb  )  ;  ■=  0  (41) 

nm  n  n  nm  n  n  nm  n  n  ■  nm  n  n 

where  kr^  is  the  argument  of  all  the  spherical  Bessel  functions  and 

Che  fact  that  h^^^  =  h^^^  has  been  used.  With  the  Wronskian  relations 


j  h«).  -  j.  b«) .  r  b<»  -i  b(»  .-i- 

n  n  -“n  n  •'n  n  ■'n  n 


h(l)  ,(2)*  .,(!)'  ,(2) 

n  n  n  n  2 

z 

where  z  is  the  argument  of  the  functions  concerned, 
we  find  that  (41)  becomes 


2i  :  :  {'c®  i2(a  +a  +2|a  '^nm  ‘ 

_  ‘  '■-nm  nm  nm  '  nm  ■  'vnm  nm 

n«0  m*0 


which  may  also  be  written  as 


n*0  m«0 


{IC®  12(1-' 2a  +1:2)  +  is®  .2(i_  2b  +1'2)}  *  o  . 
'  ^nm  •  nm  '  i-nm  nm 


It  then  follows  from  the  assumptions  on  coefficients  (33)  and  (34) 


a  C®  -  b  S®  -  0, 
nm  '■.nm  nm  •vnm 


hence  from  (39) ,  we  have  Sg  w  «  0,  0  <  r  <  min  r  but  the  analyticicy 

^  qeSD 


1 


of  solutions  of  elliptic  equations  then  implies  that  S  w  «  0, 

g 

P  £  D_u3D\{0}.  Since  -  O^p  €  3D  it  follows  then  chat 

S^w  «  0,p  c  3D  and,  with  Theorem  3.2,  Chat  k  is  an  eigenvalue  of  the 

interior  Dirichlec  problem. 

It  is  noted  that  the  sense  of  the  inequalities  in  (33)  and  (3A) 

could  be  reversed  without  changing  the  conclusion,  however  the 

sense  of  the  inequality  must  be  the  same  for  all  nonvanishing  a^ 

and  b  .  It  should  also  be  mentioned  that  inequalities  similar  to 
nm 

(33)  and  (34)  were  derived  by  Ursell  ^21  although  he  observed  that 
by  requiring  Re  a^  >  0  and  Re  b^  >  0  would  guarantee  compliance  and 
thus  attention  could  be  restricted  to  real  values  of  the  coefficients. 
These  inequalities  have  also  been  employed  recently  by  Martin  who 
explored  the  relationship  with  the  null  field  method. 


Theorem  3.3  establishes  chat  all  characteristic  values  of  -K, 

X 

are  eigenvalues  of  the  interior  Dirichlet  problem.  It  remains  to 
show  that  eigenvalues  of  the  interior  Dirichlet  problem  are  not 
necessarily  characteristic  values  of  -K^  with  g  suitably  chosen.  First 
we  establish 

Theorem  3.4:  If  g(P,Q)  is  as  in  Theorem  3.3  and  i  2a  +!!<  ,'|2b  +1  <! 

all  n,m  then  -K^  has  no  characteristic  values. 

Proof.  Assume  k  is  a  characteristic  value  so  there  exists  a  non- 
.  .  o  . 

trivial  w  which  satisfies  (35) .  The  assumptions  on  a  and  b 

nm  nm 

together  with  (45)  imply  °  all  n  >,  0  and  0  <  m  .S  n. 

With  (22)  and  (39)  and  the  vanishing  of  C®  and  S®  we  have 
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S^w 


"  ^  Sc!™  +  a_  C®_CP)}  + 

^.nin  ran 

n-0  m*0 


XUD  HID 


+  S®  ^S^(P)  +  b  S®  (P) }  «  0,  0  <  r  <  min  r  . 
^nm  nm  nm  nm  p  _  a 

qe3D  ^ 


(47) 


But  the  analyticity  of  solutions  of  elliptic  equations  then  implies 


that 


S^w  -  0,  P  £  D_C'3D\{0}.  (48) 

Hence 

rr^  S.w  -  (-1  +  Kt)w  *  0,  P  e  5D  (49) 

which,  with  (35)  implies  that  w  »  0,  violating  the  assumption,  and 
thus  establishing  the  result.  Note  that  Theorem  5.4  remains  valid 
if  the  sense  of  the  inequalities  satisfied  fay  a  and  b  is  evervwhere 
reversed. 

We  remark  that  the  three  dimensional  form  of  Ursell's  modification 
using  the  Green's  function  for  a  sphere  enclosed  by  the  scatterer 
ri3_  corresponds  to  a  particular  choice  of  the  coefficients  in  Theorem 
3.4.  It  is  however  not  necessarv  to  require  all  a  and  b  to  be 
non  zero  as  illustrated  in  the  following. 

Theorem  3.5.  If  k  is  an  eigenvalue  of  the  interior  Dirichlet  problem 

of  multiplicity  I  then  there  exists  a  g(P,Q)  with  only  4  non  zero 

coefficients  such  that  k  is  not  a  characteristic  value  of  • 

Proof.  Assume  the' contrary  so  that,  with  Theorem  3.1,  there  exists 

a  w  such  that  S.w  *  0,  pc  3D.  Observe  that  S  and  C  ^  cannot  vanish 
i  'vnm  ".•nm 

for  all  n  and  m  otherwise,  as  in  the  proof  of  Theorem  3.4,  S^w  *  0 
everywhere  causing  w  to  vanish  contrary  to  assumption.  However,  as 
shown  in  the  proof  of  theorem  3.3,  in  order  for  k  to  be  a  characteristic 
value  of  it  follows  from  (45)  that 


14. 


l?nml  ^0  1  +ll^-l)  »  0  for  all  n  and  m.  We 

choose  only  one  a^  or  to  lie  in  the  complex  plane  inside  the 
circle  of  radius  j  and  center  at  -i  when  C  or  S  is  non  zero  to 
have  a  contradiction.  If  k  is  an  eigenvalue  of  multiplicity  1  then 
the  above  argument  may  be  repeated  I  times  leading  to  at  most  i  non 
zero  coefficients.  With  such  a  g,  k  cannot  be  a  characteristic 
value  of  •  Additional  coefficients  may  also  be  taken  to  be  non 
zero  without  disturbing  the  result.  By  repetition  of  this  arginnent 
there  follows 

Theorem  3.6.  If  k^,  k2,  k^,  are  eigenvalues  (not  necessarily 

ordered)  of  the  interior  Dirichlet  problem  of  multiplicity  l.  ,  I 

N  lib 

Chen  there  exists  a  g  with  Z  i-  non  zero  coefficients  such  that 

i*l 

k^,  k,,  ...,  k^,  are  not  characteristic  values  of  -K^ ,  hence  also  not 
characteristic  values  of  and  -K^*. 

IV  The  Exterior  Dirichlet  Problem 

The  exterior  Dirichlet  problem  consists  of  finding  v_^,  a  radiating 
solution  of  the  Helmholtz  equation  in  the  sense  of  (1)  and  (2)  for 
which  v_^  »  f_^,  p  £  5D  where  f_^  is  given  in  3D.  Boundary  integral 
equations  are  obtained  from  Green’s  theorem,  (14)  and  (15)  which 
become  in  this  case 

SjW  -  (I  +  Kjt)  f_^,  p  £  5D  (50) 


and 


(I  -  K.)„.  -  ._D. 


(51) 


in  which  case  the  solution  is  represented  in  D^  via  (12)  and  (13)  as 

V  . 


1  S.w  -  i  D.f^,  P  £  D^. 

2  j  2  j  +’  + 


(52) 
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Implicit  in  this  foraulatioc  is  the  existence  of  the  normal  derivative 
of  the  double  layer.  This  can  be  assured  in  (51)  if  f_^  is  differ¬ 
entiable  on  3D  and  this  can  be  relaxed  some^diac  by  requiring  the 
boundary  to  be  Lipschitz  (Lyapunov,  of  order  1 )  [7j  .  In  the  uniqueness 
theorem  cited  below  however  there  is  always  an  implicit  assximption 
that  the  double  layer  distribution  has  a  normal  derivative  in  a 
sense  sufficient  to  allow  use  of  Green's  theorem.  That  this 
assumption  is  justified  is  proved  in  [8], 

Alternatively  the  layer  approach  starts  with  the  assumed  form  of 
the  solution 


(54) 


(55) 


v_^  *  -  DjW*,  P  £  D_^  (53) 

which  with  (10)  yields  the  boundary  integral  equation 

(I  -  K.*)  w*  =  f^.  p  £  SD. 

J 

As  before  we  denote  real  values  of  k  for  which 

(I  -  Kj)S  -  0 

has  nontrivial  solutions  as  characteristic  values  of  Kj  which  are  also 

characteristic  values  of  Kt ,  ,  and  . 

Paralleling  the  discussion  in  [3^  for  we  first  note 

.0 

Theorem  4.1.  (I-K.*)  8  »  0  if  and  onlv  if  4—  D.  w*  «  0,  o  e  aD. 

J  on  j  • 

Proof.  For  j  ”  0  the  proof  is  given  in  [Sj .  For  j  ■  1  the  argument 

0  _  o 

is  identical.  Assume  v_^  ■  -D^  w*  where  (I  -  K^)  w  ■  0.  Then  (10) 

implies  v  »  0  on  SD  and,  because  there  are  no  eigenvalues  of  the 

^  o 

exterior  Dirichlet  problem,  v_^  *  0  in  DC  3D.  Hence  4_  w*)  exists 

o  on 

and  in  fact  vanishes.  Since  w  is  an  eigenfunction  of  a  weakly 
singular  integral  equation  it  is  continuous  on  oD  and  this  fact 

together  with  the  assumption  that  3D  is  Lyapunov  of  order  1  ensures 
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3  - 

that  — ^  (Dj^  w*)  exists  and  (11)  holas  (^7]),which  establishes  the 

3n  3  2  o 

desired  result.  Conversely  if  ~  (D,  w*)  ■“  0  on  3D  let  v  =D,  w* 

3n  i  +1 

which  then  must  vanish  because  of  the  absence  of  eigenvalues  of  the 
exterior  Neumann  problem.  Since  then  vanishes  on  3D,  use  of  (10) 
completes  the  proof. 

Also,  as  shown  in  [VI  we  have 
Theorem  4.2 

-  -  3  - 

(I-K*)  w*  »  0  <=>  (D  w*)  =  0  <=>  k  is  an  eigenvalue  of  the 
o  3n  o  ° 

interior  Neumann  problem. 

As  in  Section  III,  the  fact  that  characteristic  values  of  K  are 

o 

identical  with  eigenvalues  of  tae  interior  Neumann  problem  does  not 
generalise  to  since  g(P,Q)  is  not  defined  throughout  D_.  Neverthe¬ 
less  we  can  establish  the  following. 

Theorem  4.3.  If  g  is  defined  by  (32)  subject  to  (33)  and  (34)  and 
k  is  a  characteristic  value  of  then  k  is  an  eigenvalue  of  the 
interior  Neumann  problem. 

Proof:  Assume  (I-K^)  w  =  0  has  a  nontrivial  solution  which  means  that 

o 

(I-K^)  w  ■  0  also  has  a  nontrivial  solution.  Theorem  4.1  then  implies 

,  2 

that  ■—  (D^  w*)  =  0. 

Def ine 


v  (P):*(-D^w*)(P)=-|  ^^*(q)  ^{v  (P,q)  ^  g(P,q)}ds  ,  P  £  D  \  {0} 

^  --SD  ""q  ° 

and  with  (23)  and  (32) 

o=  n  ,  • 

3  C®  (C^  (P)  +  a  C®  (P)}  + 

_  %nm  nm  nm  nm 


(56) 


V  (P) 


n*0  m»0 


+  S  ™  '•S  „(P)  ^  S  (P).',r  <  mxn  r 

^nmnm  nmnm  p  ...q 

^  qcSD  ^ 


(57) 


where  C**  and  S*  are  defined  in  (24)  and  (25)  with  ^  replacinK  w 

^TWT)  Tim  •  *  C 

Since  ■  0,  hence  ■*  0  on  D.it  follows  that,  with  Green's 

3°' 

3n- 

theorem 


i  iL  -  V  — 

J  3D  ~  3n~  "  3n 


— ^  )  ds  »  0  -  -  v_  — 

3n  ^  ^3B  ”  ‘  3n 

r 


)ds  .  (58) 

P 


Now  the  analysis  proceeds  exactly  as  in  the  proof  of  Theorem  3.3  with 
6  *  6 '  .6  6 

and  replacing  and  respectively  leading  to  the  conclusion 


e  e 

a  C  “  b  S  *  0  for  all  n  ^  0  and  0  .<  m  n. 


This  in  turn  implies  chat 


(D  w*)  (P) 

O 


r  o  a  "  n  -t  I 

i  '^<q)  g(P,q)ds„  -  E  E  (a  C®  C®  (P)  +  b  S®  S®  (P) 
I,-,  3n„  o  q  -  -  nm-vnm  nm  nm^unm  nm 

•  at>  q  n*Q  m*0 


vanishes  identically  which  means  that 


D,  w*  *  D  w*  for  all  P.  (61) 

i  o 

Sue  —  D,  w*  =  0  hence  -r—  D„  w*  »  0  and  Theorem  4.2  then  ensures  that 
3n  i  3n  o 

k  is  an  eigenvalue  of  the  interior  Neumann  problem,  thus  completing 
the  proof.  Note  that  the  theorem  remains  valid  if  the  inequalities 
in  (33)  and  (34)  are  reversed. 

While  all  characteristic  values  of  are  eigenvalues  of  the  interior 
Neumann  problem,  the  converse  is  not  necessarily  true  as  is  evident 
from  the  following: 

Theorem  4.4 

2 

If  g  is  given  by  (3.2)  and  i  a^  +  1  i  <  1,  ^  2b^+l  i<  1  for  all  n,m 
then  has  no  characteristic  values. 
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Pioof:  Assume  k  is  a  characteristic  value  so  that  there  exists  a 
o  _  0 

nontrivial  w*  such  that  (I-^)  w*  -  0.  The  fact  that  a  and  b  are 

1  nm  nm 

C  ^  6  * 

all  non  zero,  together  with  (59)  implies  that  C  •«  S  =0  for  all 

nm  nm 

n  0  and  0  4  m  n  Hence  with  (23)  and  (60)  it  is  evident  that 
o 


D,  w*.=  0,  0<r.$  min  r  . 
^  P  ,C3D  ’ 


(62) 


But  the  analyticity  of  solutions  of  elliptic  equations  then  ensures 


that  ^  ■  0,  Ft  D_\{0}  from  which  it  follows  that 


and 


(I  +  K*)  w*  =  0 


—  D,  w*  *  0. 
9n  i- 


(63) 

(64) 


But  (64)  together  with  Theorem  4.1  imply  that  (I-K*)  w*  »  0  which, 
with  (63) ,  guarantees  that  w  *  0  thus  violating  the  original  assumption 
and  establishing  the  theorem.  As  in  the  exterior  Neumann  problem 
useful  results  obtain  even  with  only  a  finite  number  of  non  zero 
coefficients  in  the  representation  of  g.  , 

Theorem  4.5. 

If  k  is  an  eigenvalue  of  the  interior  Neumann  problem  of  multiplicity 
I,  then  there  exists  a  g(P,Q)  with  only  £  non  zero  coefficients  such 
that  k  is  not  a  characteristic  value  of  . 

Proof:  Assume  that  k  is  a  characteristic  value  of  .  Proceeding  as 

6  ^  6  * 

in  the  nroof  of  Theorem  4.3  it  follows  that  a  C  =  b  =  0,  as  in 

nm'x.nm 

6  ^  0  ^ 

(59).  But  C  and  S  cannot  vanish  for  all  n>0,  0<m<nor  else 
xnm  'v.nm  ^ 

as  in  the  proof  of  Theorem  4.4,  w*  would  vanish  identically.  By 

6  *  0  ^ 

choosing  a  or  b  different  from  zero  for  one  case  when  C  or  S 

®  nm  nm  ^v.nm  -v^nm 

is  unequal  to  zero  contradicts  (59)  showing  that  k  is  not  a 
characteristic  value  of  K.  If  k  is  an  eigenvalue  of  multiplicity  ^ 


A  I'  II  ill  iilai 
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Chen  Che  argumenc  may  be  repeaced  1  cimes  vich  linearly  independenc 
w^,  w, ,  w  leading  Co  at  mosc  I  non  zero  coeff iciencs .  Repecicion 

of  this  argumenc  for  different  eigenvalues  establishes 
Theorem  4.6. 

If  . . . ,  k^  are  eigenvalues  (not  necessarily  ordered)  of  the 

interior  Neumann  problem  of  multiplicity  1,,  ...,  1  respectively 

N  l  -  N 

then  there  exists  a  g  with  I  i..  non  zero  coefficients  such  that 

i*l  ^ 

k,  ,  ...,  k^,  are  not  characteristic  values  of  ,  hence  also  not 
characteristic  values  of  K?,  and  K*. 

V  Optimal  Modifications 

Thus  far  we  have  shown  that  by  modifying  the  Green's  function  by 
adding  a  radiating  term  of  the  form  given  by  (32)  with  coefficients 
subject  only  to  the  inequalities  (33)  and  (34),  the  characteristic 
values  of  may  be  removed.  Now  we  consider  the  question  of 
choosing  the  coefficients  more  specifically  so  as  to  satisfy  various 
desirable  criteria.  One  obvious  criterion  is  to  choose  the  modified 
Green's  function  to  be  the  exact  Green's  function  for  the  problem  if  there 
exists  a  coefficient  choice  which  will  accomplish  this,  that  is,  if  the 
functions  (P)  C®  (Q)  ,S®  (P)  S®  (Q)  }  are  complete  in  D  (^  D 

Tne  task  of  finding  such  coefficients,  which  is  equivalent  to 
determining  the  Green's  fiinction  for  the  problem  is  clearly  fcrmidable. 
However  by  modifying  the  requirement  somewhat  we  arrive  at  a  coefficient 
choice  which  is  optimal  in  one  sense  of  approximating  the  Green's 
function. 

Consider  the  Dirichlec  problem  first  and  denote  the  modified 

D  C 

Green's  function  in  this  case  by  -,,v?,q).  Note  that  if  .  were  the 
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exact  Green's  function  It  would  vanish  for  all  q  t  9D  and  P  c  D  . 

^  + 

In  particular  it  would  vanish  on  any  sphere  enclosing  9D.  Thus 
Theorem  5.1:  The  choice  of  coefficients 


(C^  ,C®  ) 
nm  nm 


and  b 


'  ^^nm^  'L-(aD) 


(S^  ,S®  ) 

°°  ™  BP 

Ms®  112 

' '°nm' 'l2(3D) 


minimizes  the  ouantit'^ 


where 


(u,v)  :  *  !  u(q)v(q)ds 

J  3D 


i  i  •' 'L,(3D)  •  " 


This  choice  of  coefficients  minimizes  the  difference  between  the 

modified  Green's  function  and  the  exact  Green's  function  on  L,  (3D)  x 

L_  (SB,)  for  any  A  max  r  . 

^  ,e8D  ’ 

Proof:  Using  the  definitions  (21)  and  (32)  it  follows  that  for 


A  >■  max  r 

C 

qeoD 


! >  !  I  ^  ds 

1  L-(3D)  P 


'  /•  .  “  n 

■  ;  ;  :  :{C®  (P)  Fc^  (q)+a  C®  (q)"  + 

'  «  „  _  nm  “  nm  ^  nm  nm  ^  “ 

'  r  A ' 3D ' n*0  m*0 
p*A 


+  S®  (P)  [S^„(q)+b  S®  (q)j}|  dsds^.  (66) 

tun  nm  nm  nm  '  P 

6  0 

With  the  orthogonality  of  C  (P)  and  S  (P)  on  r  *  A  this  becomes 
®  •'  nm  nm  '  p 


n-0  nw 


-  I  S^  (q)!^}dS 

nm^  mnnm^  'nm^  nm  nm  ^  q 

.  2kA  :  r  |h^^\kA)i2{jlc^  :;,  (C®  ,cM  +r”(c"  C®  ) 

•n»0  in*0  ^  ■  ■  urn  'L^CoD)  nm  nm  nm  sD  ®nm''^nm’^nm^  3D 


■^i^nm'"'  '^nm'  '] 


L,(3D)‘^!  I^nm'  'L^(3D)‘"'^nm^^nm’^nm^  3D 


+  !b  ’  2  I  :  s®  :  ^  'i 

'  nm'  ' '  nm  'L2(3D) ' 


which  may  also  be  written 


1  D  I  2 

:  V,  ,  IT  ,  “  2kA 

,A  ^  P 

p 


:  l|h^^^kA)'2(  ,c"  '2  - IgL.-S?,.£P. 

n=0  m-0  °  4^'°^  '  C®  2 

nm'  L(,5D) 


(ci,„.c®  )  -  (S^  ,S®  )  ■ 

sn  I  ,•  nm’  nm 

a  +  ■  +  '  c-  :2  oD 

,(5D)j  nm  i:  e  .12  ,  '  nm'  L,(3D)  ,  e  2 

'  '  nm'  ‘L2(3D)  '^nm'  'L-,  (3D) 


'’'nm' 'L.  (5D)  l^nm"^  ,  ■  e 


+  i  I  S  i  :  "■  lb  + 

^  I  y^TTi  ‘  T  { 


(S^,S®  ) 
nm  nm 

oD 


nm'  ' Lo  (3D)  ;  nm  .  1  ,.e  ■  .  - 

"  ' ■  nm'  “l,(3D) 


and  this  is  clearly  minimized  when  a  and  b  are  chosen  to  satisfv 

nm  nm 

(65) ,  thus  establishing  the  theorem. 

A  natural  question  is  whether  the  choice  of  coefficients  in 
rheorem  5.1  satisfies  conditions  (33)  and  (34)  ensuring  chat  K,  has  no 
characteristic  values.  This  is  in  fact  the  case  which  we  state  as 
Theorem  5.2. 


If  k  is  real,  3D  is  not  a  sphere,  and  a  and  b  are  chosen  as 

J  r  >  jjjjj 
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in  Theorem  5.1  then  1 2a  +  1  !<  land  i  2b  +1  '<1. 

nm  nm 

Proof:  With  a  defined  by  (65)  we  have 


nm 


2(C^  ,C®  ) 
nm  nm 


(C®  ,C®  )  -2(C^  ,C®  ) 


2a  +1 
nm 


3D 


nm  nm 


1  = 


nm  nm 

3D 


|C 


nm‘ 'l,(3D) 


!  iC®  '  ' 2 

'  nm'  L2(3D; 


But  if  k  is  real, 


(69) 


C®  (q)  -  2C^  (q)  -  i  C®  (q) 
nm  ^  nm  ^  nm  ^ 


(70) 


hence 


1 2a  +1  = 

'  nm  ' 


(C®  ,C®  )  ' 

nm’  nm 

aD 


1  iC 


nm-  'L.,  ( 3D) 


(71) 


or,  denoting  by  u  and  v  the  real  and  imaginary  part  of  (-.^(q) 


; 2a  +1 , 
nm 


,  (u,u)gp>(v,v)3p-2i(u,v).p  ^ 

!  (u,u).^+(v,v).jj  I 


*  ‘  (u.u)->-(v,v))^*4  •;  (u,v)‘^-(u,u)  (v,v) } 

(u,u)  +  (v,v) 

and  this  is  clearly  less  than  one  if 

(u,v)^  <  (u,u)(v,v)  (73) 

i.e.  if  strict  inequality  holds  in  the  Schwartz  inequality.  However 
this  will  be  the  case  provided  u  and  v  are  linearly  independent  or 
explicitly  if 

=  j  (kr  )P°(cos  6  )cos  m:  (7A) 

n  q  n  a  q 

and 

=  n  (kr  )P™(cos  3  )cos  mt  (75) 

n  q  n  q  q 

are  linearly  independent  on  3D.  But  and  n^  are  linearly  independent 
solutions  of  the  spherical  Bessel  equation,  provided  r^  varies  on  3D 
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which  it  will  as  long  as  3D  is  not  a  sphere.  This  establishes  the 
desired  inequality  for  a^  and  precisely  the  same  procedure  establishes 
the  same  result  for  b^^.  When  3D  is  a  sphere  the  inequalities  (33) 
and  (34)  are  not  satisfied .However  in  that  case  explicit  results  are 
available  to  show  that  again  has  no  characteristic  values.  This  is 
demonstrated  in  Section  VI. 

A  similar  result  is  available  if  we  require  the  modified  Green's 

N 

function,  now  denoted  by  (P.q)  to  approximate  the  Green's  function 
for  the  Neumann  problem.  Specifically  we  have 
Theorem  5.3 

The  choice  of  coefficients 


(-2_  c  C  ) 

5nq  Pin»  3n  ^  nm^ 

'  i  r®  '  ' 

‘ ■  9n-  ot' 'L-(SD) 


and  b  * 
nm 


'—2—  s‘  — ^  S  ) 
■3nq  nm  in^  nm 

!  — L_  s®  '  '  2 

'  '  3nq  nm  '  L^OD) 


minimizes  the  quantity 


'1 

T —  I  ds  for  every  A  ^  max  r  . 

.nq  '  L2(3D)  p  -  q 


The  proof  of  this  theorem  is  precisely  the  same  as  the  proof  of 

Theorem  5.1  with  C®  (q)  and  S®  (q)  replaced  by  — ^  C®  (q)  and 
nm  nm  ^  ^  dnq  nm 

T -  S  (q)  and  similarly  for  C  and  S  .  The  question  again  arises 

dnq  nm  ^  nm  nm 

as  to  whether  the  coefficients  specified  in  (76)  satisfy  conditions 

(33)  and  (34)  guaranteeing  that  -K^  has  no  characteristic  values.  A 

procedure  exactly  analogous  to  the  proof  of  theorem  5.2,  using  the 

definitions  of  a  and  b  given  in  theorem  5.3  shows  that  the  inequality 
nm  nm 

i2a  +1  <  i  is  satisfied  provided  that 

'  nm 
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r-^j'  (kr  )P™(cos  6-)cos  m(f)  }  and  )P°'(cos  9  )cos  m  ) 

3nq  ■'n  q  n  q  q  3nq  n  q  n  q  q 


are  linearly  independent  on  3D.  While  this  appears  reasonable  when 
3D  is  not  a  sphere,  especially  in  view  of  the  comparable  result  in 
Theorem  5.2  involving  the  undifferentiated  functions,  a  rigorous 
proof  has  thus  far  eluded  the  authors.  Once  established,  a  similar 


procedure  would  show  that  the  inequality  i 2b^  +1  <i  is  also  satisfied. 


Yet  another  criterion  for  choosing  the  coefficients  in  the  modific¬ 
ation,  and  for  some  purposes  perhaps  the  most  meaningful  criterion  is 
to  choose  the  coefficients  so  as  to  minimize  a  bound  on  the  spectral 


radius  of  the  modified  boundary  integral  operator,  K^.  This  con' 


sists  of  minimizing 


K,  i  I  and  is  equally  applicable  in  both 

Dirichlet  and  Neumann  problems.  An  algorithm  which  accomplishes 

this  may  also  be  derived.  For  this  analysis,  however,  it  is 

convenient  to  change  notation  slightly  since  the  representation 

of  the  modification  as  as  double  sum  (32)  is  somewhat  awkward 

for  this  purpose.  Note  that  each  pair  of  integers  (n,m),0^m<n 

uniquelj'  determines  an  integer  Z  by 

2 

i  :  =  —  +  m,  0  $  m  n  (77) 


and  moreover  each  Z  >■  0  uniquely  determines  a  pair  (n,m)  by  the  same 
relation,  since  Z  lies  between  the  sum  of  the  first  n  and  the  first 
n  1  integers  for  -some  n.  With  this  relationship  between  the 


integers  Z  and  the  ordered  pairs  (n,m)  any  series  of  the  form 
oo  n  * 

A _  may  be  rewritten  as  2  A  .  Moreover  by  defining 


n*0  m*0 


nm 


£-0 


''22.1 


S*  (?) 
nm 


(78) 
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We  may  rewrite  the  modification  (32)  as 

n-0  m-0  1^0 

It  is  also  useful  for  what  follows  to  have  the  following 

,N  ^ 

Lemma  1:  The  sets  {v,.-  ^  and  {i— — }  are  linearlv  independent  on 

~  2-0  :>n  ,  -  ^ 

2-0 


3D  for  any  N  including (+  *1 

Proof.  Consider  first  the  set  {v,}' 


where  v  is  defined  in 


(78).  Assume  linear  dependence  i.e.  there  exist  {S-;  not  ail  zero 

i-0 

such  chat 


o(P)  :  -  Z  C  v^(P)  -  0,  P  e  3D. 

2*0  ^  ^ 

But  c(P),  P  £  D_^  is  a  radiating  solution  of  the  Helmholtz  equation 

with  zero  boundary  values  on  9D.  Therefore  the  uniqueness  of  solutions 

of  the  exterior  Dirichlet  problem  implies  chat  i(P)  *  0,  P  £  D  .  In 
N 

particular  *(?)  -  -  C.v.(P)  -  0  on  any  sphere  containing  9D.  But 

2-0  * 

-V  ;  are  known  to  be  linearly  independent  on  spheres  hence  C  -0^1 
violating  our  assumption  and  thus  establishing  the  linear  independence 
of  '-v  on  5d.  Similarly  the  uniqueness  of  solutions  of  the  exterior 
Neumann  problem  implies  chat  if 


■^(P)  :  *  Z  C.v,(P),  P  £  D^  and 


.  N  3v, 

1^-  :  c.  ^ 


0,  P  £  SD 


Chen  x  =  0,  P  £  D  and  again  the  linear  independence  on  spheres  implies 

3v  - 

K 

the  linear  independence  of  ;  on  3D. 


▼  ^ 
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Since  linearly  independent  on  3D,  though  not 


orthogonal,  there  does  exist  a  dual  basis  of  the  span  of  {v  , '  denoted 


.N 


by  with  the  property  that 


(81) 


In  fact  the  functions  may  be  represented  in  terms  of  {v  .;  by 

N 

^m  "  ^mi  (82) 

j«0  ■’  ■* 

where  the  coefficients  C  .  are  solutions  of  the  -equations 

mj 


(83) 


For  each  m  the  set  of  N+1  equations  is  uniquely  solvable  because  the 
linear  independence  of  {v  implies  that  the  coefficient  matrix  with 
elements  (v,,v.)  is  non  singular. 

A.  j  oU 

We  now  are  able  to  characterize  the  coefficients  of  the  modified 
Green's  function  which  minimize  the  spectral  radius  of  . 

Theorem  5.4 

The  choice  of  coefficients 


^  3D 


dD  p  r  1 
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3v , 


(84) 
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on 


I  I 


L^(3D) 


in  the  modification  (80)  minimizes  ,  ^  hence  minimizes  an 

upper  bound  on  the  spectral  radius  of 

Proof:  The  operator  norm  will  be  minimized  if  the  coefficients  in  the 


modification  minimize 


I  ;V 


! 


for  each  w  z  L-(3D).  But  in  the 


'L2(3D)  “2 


present  notation 


1'— 1 


'Vi  >4(3D)-<V-V>Sd'{ 


,f  (  f 

I  i 


N  3v 


n° 


3D'3DJ3D  p 


^  a:3^/P)v;Cq  )jw(q)w(q^). 
p  £*0  p  p  £*0  p 


ds  ds  ds 
q  qi  p 


and  we  have  a  standard  problem  of  minimizing  a  quadratic  form.  The 
necessary  conditions  for  a  minimum  are  the  vanishing  of  the  gradient 

(with  respect  to  the  coefficients)  and,  since  the  may  be  complex,  the 

derivatives  with  respect  to  the  real  and  imaginary  parts  of  3. 
will  vanish  seoaratelv.  Thus  ^ 

'  f  ^  _ 

i  j  j  (P)v  (q)  •^-2-(P,q)+  1  a  v ,  (P) v  (q, ),  - 

_  3y  N  3v.  ^ 

-  IZ®  5^  a,  •^(P)v.  (q)j  )w(q)w(q^)ds^ds^  ds^  =  0 


1-0  *  “p 


m  =  0,  1,  2,  ... 


and  therefore 


r  C  C  3v  3y  N 


1-0  "  3n  ^  '  r 
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•  aD  •  ciD  p  p  1-0  '  p 

m  >  0  .  (87) 

Since  this  is  to  be  true  for  all  w  c  L2(3D)  we  must  have 
f  3v  N  3v, 

!  ^  O;  3^(P)v-Cq)jds  -  0,  m>.  0  (88) 

•  3D  p  p  1-0  '■  p 

Forming  the  inner  product  with  v^  (q)  and  using  (81)  estblishes  that  a. 


(q)ds  ds  ds 

^  q  q^  p 


I 

i 


is  given  by  (84)  . 
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It  resnains  to  show  that  this  choice  of  coefficients  provides  a 
minimum.  That  is,  if  we  denote  by  the  modified  operator  with 
coefficients  as  specified  by  (84)  and  by  the  modified  operator  with 
any  other  choice,  we  must  verify  that 

i  |K^wj  I  <:  I  iKj^wj  I  for  all  w  £  (3D). 

Let  the  coefficients  in  the  modification  be  denoted  by 


(89) 


where  a.  are  defined  by  (84).  Then 

ilvii' 
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X  ds  ds  ds 
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(90) 


3D  3D  “  3D 

where  the  linear  terms  in  t,  vanish  because  of  the  choice  of  a° .  Upon 

K  A. 

making  the  substitution 

the  validity  of  our  result  is  seen  to  depend  on  the  positivity  of  the 


quadratic  form  j, 


3v ,  3v 
*  ,  K  IP\ 

“  ^£  ^m‘'“  ’  3n  ^ 


£■=•0  m*0 


3D* 


3v 


But,  constructing  an  orthonormal  set  from  ,  e.g. 


£-0 


V. 

using  a  Gram-Schmidt  procedure  and  the  linear  independence  of 
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(Lennna  1),  there  exists  a  set  of  coefficients  {C,.}  such  that 

-^1 
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9"  i-O  " 


for  each 


Then 
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2 
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mi 
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(91) 


where  C  is  the  matrix  with  elements  C..  and  C*  is  the  Hermitian  conjugate. 

A.* 

However  CC*  is  positive  semidef inite  jj.4,  P69'’  which  completes  the 
proof . 


ay 


We  remark  that  if  -  ,  eL^OD),  which  is  true  in  IR^  but 


9n 


not  in  IR  - ,  one  could  also  choose  the  coefficients  to  minimize 
"  ■  ^  1  •> 

(P.Q)''^  <is  ds  .  However  this  process  would  lead  to 
•  3D  •  3D  ^°p  P  ^ 

exactly  the  same  coefficients  specified  in  Theorem  5.4. 

It  should  also  be  remarked  that  while  Theorem  5.4  does  provide  an 
optimal  choice  for  the  coefficients  a.,  it  requires  the  explicit 

A. 

construction  of  .  While  this  is  possible  with  (82)  it  involves 

solving  (83)  for  the  (N+1)^  coefficients  C  .  and  this  certainlv 

mj 

represents  a  numerical  complication  for  large  K.  If  the  {v  ;  were 
orthogonal  on  3D  i.e.  when  D  is  a  sphere,  then 


:  I. 


(92) 


'L2OD) 


and  therefore  it  is  proposed  that  even  for  non-spherical  boundaries,  a 
reasonable  choice  is  provided  by 


30. 
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Although  the  consequences  of  this  choice  have  not  been  analysed,  it  is 
expected  to  be  useful  for  nmerical  purposes  at  least  for  boundaries 
which  are  not  severe  perturbations  of  a  sphere.  An  even  simpler 
approximate  choice  is  given  in  section  VII,  motivated  by  further 
results  for  a  spherical  boundary. 


The  explicit  coefficient  choices  found  in  the  last  section  simplify 
considerably  when  the  boundary  is  a  sphere  and  these  results  are 
presented  here.  They  provide  the  basis  of  a  coefficient  choice  which, 
while  not  optimal,  is  convenient  and  hopefully  useful  for  numerical 
purposes . 

Approximate  Dirichlet  Green's  Function 

It  is  perhaps  not  surprising,  since  the  modification  of  the  Green's 
function  is  in  terms  of  spherical  wave  functions,  that  in  the  case  when 
k  is  real  and  9D  is  a  sphere  the  coefficients  in  Theorem  5.1  render 
the  exact  Green's  function,  or  more  precisely  (because  of  our  choice 
of  free  space  Green's  function  (3))  differs  from  the  exact  Green's 
fxinction  for  the  Dirichlet  problem  by  a  factor  of  ■^.  If  3D  is  a  sphere 

of  radius  a  then  with  (16)  it  follows  that  for  k  real, 
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where 


cos  Sp^  -  cos  Ap  cos  5^  +  sin  5p  sin  5^  cos(Cp-:^_) 
and  the  relation 
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has  been  used.  Furthermore 
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With  (96)  we  see  that 
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therefore  using  Che  Green  theorem  approach,  the  solution  of  the 
Dirichlet  problem  is  simply,  with  (52)  and  (101) 
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On  the  other  hand  using  the  Layer  approach  we  have 
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because  the  spherical  harmonics  are  complete  on  L, (3D).  Therefore 
(5A)  becomes 

2w*  =  f_^  (104) 

which  with  (53)  again  yields  the  solution  (102) . 

Approximate  Neumann  Green’s  Function 

In  a  manner  similar  to  the  above,  when  3D  is  a  sphere  of  radius  a 
the  coefficients  specified  in  Theorem  5.3  reduce,  with  (16)  and 
orthogonality  to 
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while  for  r  ■-  a,  7^(P.q)  exists  as  a  distribution  on  L-OD) 


with  the  representation. 
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With  (107)  we  see  that 
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Therefore  using  the  Green’s  theorem  approach  the  solution  of  the 
Neumann  problem  is  with  (28)  and  (109) 
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On  the  other  hand  using  the  layer  approach  we  have,  with  (108)  and 
the  completeness  of  the  spherical  harmonics  on  L^C  D) 
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Therefore  (30)  becomes 
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which  with  (29)  again  yields  the  solution  (110) 


(112) 


Minimum  Snectral  Radius 


It  is  of  some  interest  to  compare  the  above  results  with  the 
coefficients  obtained  by  minimizing  the  spectral  radius  of  the 
modified  operator  when  5D  is  a  sphere.  The  coefficients  specified 
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in  Theorem  5.4  simplify  as  noted  in  (92)  and  (93).  Thus  with  (78), 
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A  similar  analysis  shows  that  in  this  case 
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With  a  and  b  so  defined  the  modified  Green's  function  becomes,  using 
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(95)  and  (106) 
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Hence  with  (98)  (99) ,  (106)  and  (107)  we  see  chat 
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Moreover  with  (100)  and  (108)  we  have 
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But  this  last  result  implies  that 

K-j_w  »  ic,^  -  0  ^  w  £  L,(3D)  (119) 

hence  the  spectral  radius  of  the  modified  operator  is  zero.  Observe 
that  with  (103)  and  (111)  the  spectral  radii  of  the  modified  operators 
generated  by  vj  and  are  both  equal  to  one. 

As  a  consequence  of  (11^)  we  have 

Dl(Sig^)  *  -  S^g^,  P  e  £  L2(3D)  (120) 

and 
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Equation  (120)  follows  from  the  uniqueness  of  solutions  of  the  exterior 
Dirichlet  problem  and  observing  that  with  (10)  and  (119) 
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whereas  (121)  follows  from  uniqueness  for  the  exterior  Neumann  problem 
and  observing  that  with  (9)  and  (119) 
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To  solve  the  Dirichlet  problem  with  this  modified  operator  we  have 
with  (51)  and  (119)  and  (52) 
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but  with  (121)  and  (117)  this  becomes 
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in  agreement  with  (102).  Similarly  using  the  layer  ansatz, 
(119)  imply  that 


(125) 
(54)  and 
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and  then  (53)  yields  the  same  solution  (125) . 

To  solve  the  Neumann  problem,  (26)  and  (119)  yields 
w*  -  S^g_^ 


and  with  (28) 
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Using  (120  and  (116)  this  becomes 
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in  agreement  with  (110).  With  the  layer  approach  (30)  and  (119)  imply 

w  =  g^  (129) 


which  together  with  (29)  again  yields  the  solution  (128)  . 

The  remarkable  result  (119)  shows  that  the  spectral  radius  of  the 
modified  operator  is  zero  in  this  sphere  example  when  none  of  the 
coefficients  in  the  modification  vanish.  For  applications  one  would 
like  to  know  the  effect  on  the  spectral  radius  of  the  modified  operator 
■if  only  a  finite  number  of  coefficients  are  chosen  optimally  while  the 
remainder  are  taken  to  be  zero  i.e.  modifications  of  the  form 
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coefficients  are  again  given  by  (113)  and  (111) .  However,  instead  of 
(118)  we  find  that  for  r^  “  a. 
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so  that,  using  the  orthogonality  of  the  spherical  harmonics, 
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Using  very  crude  estimates  for  the  spherical  Bessel  and  Hankel  functions 
it  may  be  shown  that  there  exists  a  constant  C,  independent  of  ka  and  N, 
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Thus  for  any  value  of  ka  it  is  clear  that  i |K^| |  <1  for  N  large  enough, 
and  therefore  the  spectral  radius  of  the  modified  operator  can  be  made 
less  than  one  with  only  a  finite  number  of  nonzero  coefficients  in  the 
modification. 

VII  Concluding  Rgrtmrks 

In  this  paper  we  have  shown  how  the  modified  Green's  function 
approach  of  Jones  and  Ursell  can  be  extended  to  both  Dirichlet  and 
Neumann  problems  for  the  Helmholtz  equation  in  three  dimensions  giving 
rise  to  an  integral  equation  of  the  second  kind  that  is  uniquely 
solvable  for  all  real  values  of  k.  In  addition  we  have  shown  how  to 
choose  the  coefficients  in  the  modification  optimally,  either  to  best 
approximate  the  Dirichlet  or  Neumann  Green's  function  or  to  minimize 
the  norm  of  the  modified  operator.  The  optimal  results  were 
exhibited  explicitly  for  the  sphere  where  it  was  also  shown  that  only  a 
finite  number  of  coefficients  need  be  chosen  different  from  zero  to 
force  the  spectral  radius  of  the  integral  operator  to  be  less  than  one 
for  any  finite  value  of  ka. 

The  coefficients  for  which  the  modified  Green's  function  best 
approximates  the  exact  Green's  function  are  given  explicitly  (Theorems 
5.1  and  5.3).  However  the  coefficients  which  minimize  the  spectral 
radius  (Theorem  5.4)  require  the  construction  of  a  dual  basis  for 
•C®  ,  S®  }  on  SD.  For  nonspherical  SD  and  large  values  of  N  this 
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may  be  a  considerable  numerical  complication  and  one  possible 
simplification  has  been  proposed,  (93) .  The  explicit  results  for 
the  sphere  allow  still  another  even  simpler  coefficient  choice. 
Comparing  the  explicit  results  for  the  coefficients  which  give  the 
best  approximation  to  the  Dirichlet  and  Neumann  Green's  functions 
when  3D  is  a  sphere,  equation  (94)  and  (105)  with  the  coefficients 
which  minimize  the  norm  of  the  modified  operator  (113)  ,  we 
find  that  these  latter  coefficients  may  be  written  as 
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While  the  coefficients  thus  defined  minimize  the  norm  of 
It,  only  when  3D  is  a  sphere,  it  is  proposed  that  this  choice  be  used 
for  nonspherical  surfaces  as  well.  It  is  easy  to  demonstrate  that 
conditions  (33)  and  (34)  which  ensure  that  has  no  characteristic 
values  is  fulfilled  since 
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where  the  last  Inequality  follows  from  Theorem  5.2.  A  similar  analysis 
shows  that  \  2b^*l  [<1.  It  is  reasonable  to  expect  that  with  the 
coefficients  defined  by  C137)  and  (138)  the  spectral  radius  of  the 
modified  operator,  which  is  zero  when  3D  is  a  sphere,  will  remain  less 
than  one  for  a  class  of  nonspherical  surfaces  and  furthermore  if  only 
a  finite  number  of  coefficients  are  taken  to  nonzero,  defined  by 
(137)  and  (138) ,  the  spectral  radius  will  be  less  one  in  an  interval  in 
k  which  increases  with  the  number  of  nonzero  coefficients.  A 
characterization  of  this  class  of  surfaces  and  relationship  between 
the  geometry,  values  of  k  for  which  the  spectral  radius  is  less  than 
one,  and  the  number  of  nonzero  coefficients  constitute  a  class  of 
problems  that  remain  to  be  solved. 
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